Conserving approximations in time-dependent quantum transport: Initial correlations 

and memory effects 
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We study time-dependent quantum transport in a correlated model system by means of time- 
propagation of the Kadanoff-Baym equations for the nonequilibrium many-body Green function. 
We consider an initially contacted equilibrium system of a correlated central region coupled to 
tight-binding leads. Subsequently a time-dependent bias is switched on after which we follow in 
detail the time-evolution of the system. Important features of the Kadanoff-Baym approach are 1) 
the possibility of studying the ultrafast dynamics of transients and other time-dependent regimes 
and 2) the inclusion of exchange and correlation effects in a conserving approximation scheme. We 
find that initial correlation and memory terms due to many-body interactions have a large effect 
on the transient currents. Furthermore the value of the steady state current is found to be strongly 
dependent on the approximation used to treat the electronic interactions. 
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The ultimate goal of molecular electronics [l[ in solid 
state circuitry is to miniaturize the size and maximize 
the speed of integrated devices. Advances in this field 
crucially depend on the accumulated experimental and 
theoretical knowledge. For the latter to progress it is es- 
sential to develop quantum mechanical approaches that 
are able to deal with open and interacting systems in an 
out of steady-state regime. Desirable features of such ap- 
proaches are therefore 1) the possibility to study the ul- 
trafast dynamics of transients and other time-dependent 
(TD) regimes and 2) the inclusion of exchange and corre- 
lation effects in a conserving approximation scheme. Fea- 
ture 1) was incorporated in some recently proposed one- 
particle frameworks and was exploited to address sev- 
eral issues in TD quantum transport (QT) 0, S, S @. 
These frameworks can, in principle, be combined with 
TD density functional theory 0, @, @, , thus providing 
a route to include Coulomb interactions (possibly in a 
conserving way [Io||). Feature 2) is an essential require- 
ment as realistic time evolutions must preserve basic con- 
servation laws as, for instance, the continuity equation. 
Conserving approximations like, e.g. self-consistent 
Hartree-Fock (HF), second Born (2B) or GW, have re- 
cently been employed in the context of QT but the im- 
plementations have been, sofar, restricted to steady-state 
regimes [H, O Q E [3 . 

In this Letter we propose an alternative approach to TD- 
QT that encompasses both feature 1) and 2). It is based 
on the real-time propagation of the embedded Kadanoff- 
Baym (KB) equations 17, 18, Qji H3] which are equa- 



tions of motion for the nonequilibrium Green function 
from which basic properties of the system can be calcu- 
lated. We consider a set {a} of noninteracting electronic 
reservoirs connected via a tunneling Hamiltonian to an 



interacting many-body quantum system C. The Green 
function G(z, z') (we suppress basis indices) projected on 
C obeys the equation of motion [3] 



id. 



■h(z)]G(z,z') = S(z,z') + J^dz T,(z,z)G(z,z') (1) 



where z and z' are time-coordinates on the Keldysh con- 
tour c 19]. We consider systems initially (times t < 0) 
contacted and in equilibrium at inverse temperature /? 
and chemical potential /i. The corresponding contour 
is described in Refs. and IbSl. In Eq.{T]), h(z) is the 
one-body Hamiltonian of the interacting system C and 
£ is the time-nonlocal self-energy. The latter describes 
the effects of many-body interactions and embedding of 
the system and is the sum of a many-body self-energy 
Smb[G] and an embedding self-energy £ C mb- The former 
is a functional of the projected Green function G only and 
can be expressed in terms of Feynman diagrams while the 
latter is a sum S em b = J2 a ^™b,a where 



J cmb. 



,{z, z') = t Ca (z)g a (z, z')t a c(z'). 



(2) 



In Eq.([2]) g a is the Green function of the uncontacted 
lead a and matrices tc a and t a c describe the couplings of 
system C to the leads. The TD equations obtained from 
Eq. (fT]) by restricting time-arguments to different parts of 
the Keldysh contour are known as the Kadanoff-Baym 
equations [U El, H, 13] and are the main tools of this 
work. As the system is driven out of equilibrium by a TD 
bias voltage, the current flowing into lead a is obtained by 
taking the time derivative of the total number of particles 
in a 2JJ and reads 

I Q (t) = -2ReTr c [G< • ££ emb + G R • S< omb ](t, t) 

-2ReTr c [GUs^ emb ](t,t) (3) 
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where the trace is taken over the central region indices 
and the products • and * denote integrations over the 
real and imaginary tracks of the contour (see Ref0 for 
details). The objects superindexed with ^, ] [ correspond 
to time arg uments on different parts of the Keldysh con- 
tour 0, UK and R/A denote the retarded and advanced 
components. The last term in Eq.Q arises from integra- 
tion along the imaginary branch (0, —i/3) of the Keldysh 
contour [3, [T3| and explicitly accounts for the effects of 
initial correlations and initial-state dependence. If one 
assumes that both dependencies are washed out in the 
long-time limit (t — * oc) then the last term in Eq.Q 
vanishes and the Meir-Wingreen formula [2l| is recov- 
ered. The KB approach provides a natural tool to inves- 
tigate the validity of this assumption which has remained 
unexplored sofar, and that we partly address below. 
Using the KB equations we first solve the embedded and 
correlated equilibrium problem and then propagate the 
system in time after applying a time-dependent bias (cf. 
Fig.l of Refll8l). For this we extend the implementation 
of Refllato open systems, i.e., by replacing Smb[G] with 
Smb[G] + S C mb- Different TD perturbations allow us to 
address several open issues in correlated TD-QT. 1) We 
can set the tunneling Hamiltonian to zero at the initial 
time or not (which corresponds to C initially contacted 
or uncontacted) and study the effects of the initial condi- 
tions (partition-free Q vs. contacting approach [23|) on 
the TD density and total current. 2) Many-body interac- 
tions can be included in the initial solution of the equilib- 
rium problem or switched on at later times. In this way 
effects of initial correlations [T3| on transient and steady- 
state properties can be highlighted. 3) The dependence 
on the history of the applied bias can be investigated for 
different approximate £mb[G]'s. Due to the nonlinearity 
of the problem and possibly to the nonlocality in time of 
the many-body self-energy nontrivial memory effects may 
occur (bistability, hysteresis phenomena, etc.). 4) We can 
study ac biases, pulses or other kind of TD biases as well 
as TD gate voltages and TD contacts in correlated QT. 
We stress that for a given approximate Smb [G] ah kinds 
of TD perturbations within the KB approach require the 
same computational effort. 

To study several of the issues mentioned in points 1)- 
4) above we consider an interacting device coupled to 
noninteracting one-dimensional leads. The full system is 
described by: 



a) 0.12 



H(t) 



ij,(TO> 



ij.crcr' 



a c jcra + c jaa 



f I. 



(4) 



where i,j are site indices, a, a' spin indices and where 
$ , c and d' , d are the creation and annihilation opera- 
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FIG. 1: (color online) Transient currents for the HF (red),2B 
(green) and GW (blue) approximations V = 0.5 (a) and V = 
0.8 (b). The spectral functions are displayed in c) (V = 0.5) 
and in the inset of b) (V = 0.8). The inset in c) depicts 
the modulus of the Fourier transform of the transients, d) 
density occupation number n(t) on site 1. The inset shows 
the modulus |n(tj)| of its Fourier transform 



tors for leads and central region respectively. The first 
term in Eq.((3]) describes the leads with a TD bias U a (t) 
while the second and third term describe the one-body 
and many-body interactions of the central region. Fi- 
nally the last term describes the coupling between the 
leads and the central region. Since we consider semi- 
infinite leads, the energy band of the leads is continuous 
with a finite band width and we therefore do not employ 
the often used wide band limit approximation. We stress 
that the KB approach is not limited to I-D leads as the 
leads enter only via the embedding self-energy. 
We consider a system with two central sites coupled to 
leads located on the left and right of the central region 
(i.e. a = L,R). We use the parameters in = t22 = 
and ti2 — t2i = — 1. For the many-body interactions 
we take vn = V22 — 2 and V12 = V21 = I. The chemi- 
cal potential fi for the whole initial equilibrium system is 
set at the middle of the HF gap and is determined by a 
Hartree-Fock calculation on the uncontacted but corre- 
lated central region which yields the value /i = 2 with the 
parameters described above. For the semi-infinite leads 
we use tf/ R = and tfj = — 1.5 if i and j are neighbor- 
ing sites and zero otherwise. The leads are coupled to the 
central region by coupling elements V\^l — V2JR = V 
when j is the first site on the left or right lead and zero 
otherwise (we use V — 0.5 and V = 0.8). We further 
take [3 = 90 (zero temperature limit). All quantities are 
expressed in atomic units. For Emb we employ the HF, 
2B and GW approximations which have also been used 
in earlier transport studies [H, [l(|. The GW approxi- 
mation includes the dynamical screening of the electron 
interaction whereas 2B includes all Feynman diagrams to 
second order in the bare interaction. In the small system 
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that we study here the second order exchange diagram 
(incorporated in 2B but not in GW) is not negligible and 
therefore the 2B approximation gives probably the most 
accurate description of electronic correlations. 
Correlations in transients. In Figfj] we show the tran- 
sient currents flowing into the right lead and the spec- 
tral functions for the HF, 2B and GW approximations. 
The system is driven by a symmetrically switched bias 
U L (t) = -U R {t) = U6{t) with U = 1.0, where 6(t) is 
a Heaviside function, i.e. we consider a sudden switch 
on of the bias at t = 0. Results are displayed for weak 
(V = 0.5) and strong coupling (V = 0.8) of the central re- 
gion to the leads. We first consider the spectral functions 
which are defined as A(T,u>) = -Tr Im / dte iut [G> - 
G<](T+|,T-|) where T= (t 1 +t 2 )/2 and t = t 1 -t 2 . 
We find that after the steady state has been reached the 
spectral functions do not depend on T anymore. In equi- 
librium the spectral peaks are at e° 2 = 0.5,3.5 for HF, 
2B and GW. For the biased system (Fig.lc and inset 
of Fig. lb) the electron correlations beyond HF lead to 
a narrowing of the gap between the spectral peaks and 
a broadening of the spectral function. This bias depen- 
dent gap closing mechanism was recently identified by 
Thygesen [l6j |. The positions £1,2 of the spectral peaks 
strongly affect the final steady state currents as they are 
largest when both spectral peaks enter the bias window, 
i.e. for biases such that fj, ± U rs ei j2 . This condition is 
much better satisfied for 2B and GW than for HF and 
explains the higher values in 2B (the highest) and GW 
(top panels of FigfT]). Let us now focus on the temporal 
structure of the transients. The transient currents show 
an oscillation that becomes more pronounced when we 
weaken the coupling from V = 0.8 to V — 0.5. The 
modulus of the Fourier transform of the current (minus 
its steady state value) is displayed in the inset of Fig.lc. 
There is a frequency peak at 2.5 in all many-body ap- 
proximations and for HF also one at 0.5 (for 2B and GW 
there is a broad peak around zero). These frequencies 
cannot be directly related to the spectral functions of 
FigfT] as those correspond to the steady state limit when 
the transients have settled. The frequencies instead cor- 
respond to transitions between the spectral peaks e? 2 °f 
the initial equilibrium system (which for HF,2B and GW 
have similar values of t\ 2 — fi « ±1.5) and the incom- 
ing/outgoing states at the left/right Fermi energy /j,±U. 
The peak at 0.5 is not visible for 2B and GW because the 
corresponding oscillation is damped faster than in HF. In 
the case V = 0.8 the current oscillations are suppressed 
and the steady state is obtained earlier compared to the 
V = 0.5 case since the electrons can tunnel in and out of 
the device more easily. The oscillations have a clear re- 
lation to density changes in the central region. In Fig. Id 
we display the TD site occupation ni(t) = —iG^it^t) 
of site 1 in the central region (the occupations on site 1 
and 2 satisfy n\{t) + n 2 {t) ~ 2). The sudden switch-on 
of the bias generates a density oscillation in the central 



region which damps on a time-scale comparable to the 
time to reach the steady state current. In this limit the 
system becomes polarized and part the electron density 
is accumulated to the right side of the device i.e. in the 
direction in which the current flows. The modulus of 
Fourier transform of n\{t) (inset) displays peaks at ex- 
actly the same frequencies as obtained from the wiggles 
in the transient currents. 

Conservation of charge. Since we use conserving ap- 
proximations particle number must be conserved in the 
system. This is illustrated in Fig. 2a for the 2B ap- 
proximation: the system is driven out of equilibrium 
by an asymmetric steplike bias Ul = 0.9, Ur = —0.4. 
This plot shows the currents Il/r{£) as well as the time 
derivative of the number of particles N(t) in the cen- 
tral region: clearly the particle number conservation law 
II + Ir = —alN/dt is obeyed. 

Initial state dependence and memory. Initial correlations 
manifest themselves in two ways in the KB equations. 
First, the initial values of the time-dependent Green func- 
tions are determined by the equilibrium Green function 
at t = that is obtained by considering both time ar- 
guments on the vertical track of the Keldysh contour. 
Second, the KB equations contain terms that describe 
memory of the initial state during the time propagation. 
These terms depend on the self-energies 
with mixed real and imaginary time arguments. We in- 
vestigate these two memory effects separately by either 
setting the self-energy Kmb(z, z') to zero for z and z' on 
the vertical track of the contour (initial state is noncor- 
related) or by setting s]^[ b and/or to zero (initial 
state is correlated but memory effects due to embedding 
and/or electron correlations are neglected). 
We start by setting Smb to zero when both time- 
arguments are on the vertical track of the contour for 
a situation in which we propagate an unbiased system 
from time t = to a finite time to at which time 
we switch on a sudden symmetric bias, i.e. we use 
Uh{t) = —Uu(t) = 9(t — t )U. Since electron correlations 
are taken into account in the time propagation but not in 
the initial state there will be a charge redistribution for 
times t > 0. The result is compared to an initially cor- 
related KB propagation. We take to = 20 and U = 1.0. 
The noninteracting system at t = has spectral peaks 
at energies e? 2 = ±1.0 (i.e. below the chemical poten- 
tial at n = 2). As a result of electron interactions for 
t > we find upward shifts in the spectral peaks yield- 
ing one peak above and one peak below fi. As a conse- 
quence a charge of about 2 electrons is pushed into the 
leads. The corresponding current, shown in Fig. 2b for the 
2B approximation, is saturated before the bias voltage is 
switched on at time to. For later times t > to the tran- 
sient currents with inclusion and with neglect of initial 
correlations are indistinguisable. We therefore conclude 
that the initially uncorrelated system has relaxed to a 
correlated state when the bias is switched on. 
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FIG. 2: (color online) a) Transient currents for asymmet- 
ric bias (see text), b) Transient currents for 2B with and 
without the initially interacting ground state (2B and 2B0 
correspondingly)^) Transient currents in HF and 2B approx- 
imations with and without the memory terms S^b/MB ( see 
text), d) Transient currents for different bias-switchings. All 
panels correspond to V = 0.5 



To study how initial states are remembered during time- 
propagation we compare full solutions of the KB equa- 
tions to ones in which we neglect the terms E^j, and/or 

Syjj. However, at the initial time t — we still employ 
the fully correlated embedded equilibrium Green func- 
tion. The results are displayed in Fig. 2c for the HF and 
2B approximations. We find that neglect of the mem- 
ory terms E^ / ^ has a considerable effect on the transient 
currents. In the HF case these terms only contain the 
embedding self-energy E^[ b ( as ^mb of HF is purely lo- 
cal in time) and therefore the term describes memory of 
the initial contacting of leads. Neglect of this term leads 
to the curve labeled HF 1 in Fig. 2c. For the 2B case 
there is also a dependency on the many-body self-energy 
Syg. We therefore have two curves for 2B, one in which 
we neg lect only E^|, (labeled 2B 1) and one in which we 



neglect both E]^[ b and (labeled 2B 2). We see that 
neglect of E^[ b has a considerable effect on the transients 
while neglect of only E^g has a smaller but still notica- 
ble effect. We further see that the same steady state 
current develops as with the memory terms included and 
therefore conclude that the memory terms eventually die 
out in the long-time limit. This is in agreement with the 
memory loss theorem proven in RefsQ andd for the case 
of Green functions that are sufficiently smooth. We fi- 
nally note that there are situations for which the Green 
function is not a smooth function in which case persis- 



tent oscillations may appear, see Refl24 . 
Time dependence of applied bias. We finally investi- 
gate the dependence of the transient currents on various 
forms of the time-dependent bias. In Fig. 2d we show the 
2B transient currents driven by different TD symmet- 



ric biases: U L (t) = -U R {t). We take U L (t) = U0(t), 
U L (t) = [/Erf(wii) and U L (t) = C/sin 2 (w 2 i) for t < ^ 
and U L (t) = U for t > ^ with U = 1.0, lu 1 = 0.5 
and u>2 = 0.1. We observe that the sudden switch-on 
produces rapid oscillations. They are more damped with 
slower switch on of the bias voltage. The steady state 
currents are, however, the same for all three cases. How- 
ever, due to nonlinearity of the KB equations existence of 
bistable solutions and hence different steady states may 
be possible for different biases. This will be part of future 
investigations. 

We conclude that the KB equations provide a powerful 
tool to study correlated quantum transport in real time. 
The method allows for inclusion of many-body correla- 
tions while satisfying important conservation laws. We 
found that many-body interactions have large effects on 
steady-state and transient currents. The temporal fea- 
tures in the transients and density distributions were an- 
alyzed in detail and related to level structure displayed 
in the spectral functions, a study of utmost importance 
to interpret transport spectroscopy experiments. We fur- 
ther showed that memory terms have large effects on the 
time-dependent currents. 
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